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ABSTRACT 

We have analyzed high spatial, moderate spectral resolution observations of Eta Carinae (77 Car) 
obtained with the Space Telescope Imaging Spectrograph (STIS) from 1998.0 to 2004.3. The data 
were obtained at discrete times covering an entire 2024 day spectroscopic cycle, with focus on the X- 
ray/ionization low-state which began in 2003 June. The spectra show prominent P-Cygni lines in H I, 
Fe II and He I which are complicated by blends and contamination by nebular emission and absorption 
along the line-of-sight toward the observer. All lines show phase and species dependent variations in 
emission and absorption. For most of the cycle the He I emission is blueshifted relative to the H I and 
Fe II P-Cygni emission lines, which are approximately centered at system velocity. The blueshifted 
He I absorption varies in intensity and velocity throughout the 2024 day period. We construct radial 
velocity curves for the absorption component of the He I and H I lines. The He I absorption shows 
significant radial velocity variations throughout the cycle, with a rapid change of over 200 km s _1 
near the 2003.5 event. The H I velocity curve is similar to that of the He I absorption, though offset 
in phase and reduced in amplitude. We interpret the complex line profile variations in He I, H I and 
Fe II to be a consequence of the dynamic interaction of the dense wind of 77 Car A with the less dense, 
faster wind plus the radiation field of a hot companion star, 77 Car B. The companion's wind carves 
out a cavity in 77 Car A's wind, which allows UV flux from 77 CarB to penetrate and photoionize an 
extended region of 77 Car A's wind. During most of the orbit, 77 CarB and the He + recombination zone 
are on the near side of rj Car A, producing blueshifted He I emission. He I absorption is formed in the 
part of the He + zone that intersects the line-of-sight toward r\ Car. We use the variations seen in Ho I 
and the other P-Cygni lines to constrain the geometry of the orbit and the character of 77 Car B. 

Subject headings: stars: binaries:spectroscopic, stars: individual: Eta Carinae, stars:winds 



1. INTRODUCTION 

Eta Carinae (77 Car) is a super luminous, unstable 
object which underwent a rapid brightening in the 
1840s, accompanied by ejection of a substantial amount 
of material. These ejecta now form a structured, 
bipolar nebula, the Homunculus, which makes direct 
observations of the star difficult. Hence, the nature of 
77 Car has long been debated. Most observations up 
through the mid 1990s could be interpreted as arising 
from either a single star or binary system (| Davidson! 
Il997| ) . More recent evidence favors the binary star inter- 
pretation. Indirect signatures of the companion star are 
periodic events first discovered as intensity variations 
in He I A10830 (|Daminelilll996l . and references therein) 
with a 5.5 year periodicity. The excitation events are 
accompanied by eclipse-l ike minima observed with U BV 
and BVR photometry (Ivan Genderen et all |2003|). in 
near infrared JHKL phot ometry (Whitcl ock et al.ll2004h 
and in X-ray brightness (jlshibashi et al.l["l999l : ICorcoranl 
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2005!) . TheX -ray light-curve as measured with the Rossi 
X-ray Timing Explorer (RXTE) is characterized by a 
gradual increase in brightness just before a rapid decline 
to a low-state which lasts ~70 days. RXTE observations 
of two such cycles fixed the orbital period to 2024 ± 2 
days (5.54 years), which is in agreement with the period 
derived from the infrared observations and from the 
He I A10830 variability. Emissio n from the Weigelt 
blobs (jWeigelt fe Ebersberger1ll986| ) located within 0'.'3 
of 77 Car also shows low-excitation intervals tracing the 
spectroscopic period, with emission from highly ionized 
species such as [Ar III] and [Fe IV] disappearing during 
the short spectroscopic minimum. IVerner et alJ (2005) 
attributed this variation to the visible presence of an 
additional source of flux, possibly an O or Of/WN7 
star with T e ff~37,000 K, during the time outside of the 
minima. 

The X-ray variations, as seen by RXTE, Chandra and 
XMM-Newton, have been modeled as a collision between 
a massive wind of the primary star (77 Car A) and a less 
dense wind of th e hot companio n (77 C a rB) in a highly 
ellipt i cal orbit jlshibashi et al.1 119991 : ICorcoran et al.1 
2001} lHamaguchil I2007I l Additional evidence for the 
presence of a binary component is the detection of He II 
A4686 which appeared in g round-based spectra of 77 Car 
(jSteiner fe Daminelil l2004t IStahl et all l2005h shortly 
before the mini mum but d isappe ared at the o n set o f 
the minimum. iGullI (|2005l ) and iMartin et all (|2006f ) 
confirmed the presence of He II A4686 with high angular 
resolution Hubble Space Telescope (HST) STIS spectra, 
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and lMartin et all (|2006f ) proposed that the He II A4686 
emission is formed in the interface between the stellar 
winds or possibly arising from photo-ionization by the 
companion star. 

The multi-wavelength variations observed in 77 Car 
have been interpreted as arising from the interaction 
of the winds of 77 Car A and 77 CarB or the impact of 
the radiation from the hot 77 Car B o n the circumstellar 
material. Recently, llping et ail (|2005l ) reported evidence 
of the companion through detection with the Far 
Ultraviolet Spectroscopic Explorer (FUSE) of a source 
of far-UV radiation which disappeared just before the 
beginning of the 2003.5 event, and reappeared by 2004 
March. 

Binary syste m mode ls have been dis cuss ed b y 
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(119971): llshibashi et al 



( 20011): IFalceta-Goncalves et all 



1999); 



2005) 



Soked (20051) ? among others. In these models the far-UV 



flux from a hot companion alters the ionization state 
of 77 CarA's wind, while X-rays are produced in a 
bow-shock due to the collision of the dense wind of 
77 Car A with the lower density, higher velocity wind of 
r\ CarB. In most of the colliding wind binary models, 
the low-state in observed X-ray intensity, the brightness 
in He I A10830 and in the visible-band, IR and radio 
fluxes, occurs when 77 CarB along with the wind-wind 
collision region moves behind 77 Car A and is occulted 
by its dense wind. Other interpretations are that the 
reduction in X-ray flux is primarily due to mass transfer 
between the stars for a brief period near periastron 
passage (jSoker 2005) or that the X-ray flux is reduced by 
a dense accumulation of wind material trailing 77 CarB 
after periastron passage (j Falceta- Gonpalves et all 



2005). In the Falc eta-Goncalves et al 
companion star is positioned between 



model, the 
the observer 



and the primary star during the periastron passage. 
This orientation is difficult to reconcile with the 
radiation cut-off observed toward the Weigelt blobs 
during the minimum. The Weigelt blobs are located 
on the same side of 77 Car as the observer, based 
on pr oper motions and obse r ved bl ueshifted emis- 
sion (iDavidson fc Humphreys! Il997t IZethsonl 120011 : 
iNielsen et al.ll2007l) . The IFalceta-Goncalves et al.l orien- 
tation is also difficult to reconcile with HST Advance d 
Camera for Surveys (ACS) imagery (|Smith et al.l l2004). 

The 77 Car binary system bears striking similarities 
to the canonical long-period colliding wind binary, 
WR 140. WR 140 is a massive system (total system 
mass ~70M Q ) composed of a WC7 +04-5 pair in a 
7.9-year, highly eccentric orbit (Marchen ko et al1l2003f) . 
The strong wind from the WC7 star dominates the 
emission from the system, though the 04-5 star is 
the more massive component and dominates the total 
luminosity of the system. WR 140 is an X-ray variable, 
with an X-ray light-curve which rapidly climbs to a 
maximum of ^5xl0 33 ergs s _1 in the 0.5—10 keV 
band, followed by a rapid decline to minimum near 
periastron passage ([Pollock et al.l 120051) . L ike 77 Car, 
WR 140 shows infrared ( Williams et al.l fl990) and radio 
(|White fc Beckerl I1995D variability correlated with the 
orbital period. Like 77 Car, the infrared brightness of 
WR 140 peaks near periastron passage, while the radio 
flux drops precipitously at that time. Unlike 77 Car, 
in WR 140 both stellar components can be directly 




triplet 



singlet 



Fig. 1. — Energy level diagram including all He I transitions used 
in this analysis. Note, the amount of energy needed to populate 
the excited states (1 cm -1 = 1.240xl0 -4 eV). All wavelengths are 
in vacuum. 



observed; and unlike 77 Car, WR 140 is not shrouded by 
thick circumstellar ejecta. 

As part of an HST Treasury project 7 , 7; Car was 
observed with HST STIS throughout the 2024 day (5.54 
year) cycle including the spectroscopic low-state which 
began on 2003 June 29. The STIS spectra provide the 
highest spatial resolution spectrometry of the star ever 
obtained and thus offer the best data-set to decouple 
stellar changes from changes in the circumstellar nebula. 
STIS's extended wavelength coverage into the ultraviolet 
probed for direct evidence of the companion star in 
the form of excess continuum and potential wind lines 
characteristic of a hot star. The optical wind lines in 
the stellar spectrum are excellent diagnostics of the 
companion's influence on 77 CarA's wind. Of particular 
importance are the P-Cygni lines from excited levels in 
neutral helium. The population of the excited helium 
energy states requires high energy photons, and is 
dependent o n the radiation from the companion star. 
IVerner et al.l (|2005l ) found that radiation from 77 CarB 
is necessary to produce nebular He I emission lines 
observed in the spectra of the Weigelt blobs. 

This paper describes the characteristics of the He I 
lines in the STIS spectra throughout the orbit, and 
compares their behavior to lines in H I, Fe II and [N II]. 
Since these lines are formed under different physical 
conditions in the 77 Car system, they sample the response 
of the system to the interactions with the companion.Wc 
derive a velocity curve based on the He I absorption, 
and show that it is similar to the expected photospheric 
radial velocity curve for a star in a highly eccentric 
orbit. The spectral line analysis is used to constrain 
the system parameters and indicates that the He I 
absorption column samples a spatially complex region 
that changes continuously throughout the spectroscopic 
period due to dynamical and radiative interactions 
between the two stars in the 77 Car system. 

7 see http://archive.stsci.edu/prepds/etacar/ 
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TABLE 1 

Data used in this analysis. Observations made with the G430M or G750M grating. 



Proposal id. 


Observation Date 


JD 


Orbital Phase, <f 


) a Position Angle 






(+2,450,000) 




(dog) 


7302 


1998 January 1 


0814 


0.007 


— 100 




199o Marcn ly 


0891 


0.045 


—28 


QflQfiC 
oUoD 


lyyo iNoveniuer zo 


1 1 A O 

114z 


U. lOy 


i n 
loo 




1999 lebruary 21 


1231 


0.213 


—28 


8327 b 


2000 March 13 


1616 


0.403 


-41 


8483 c 


2000 March 20 


1623 


0.407 


-28 


8619 b 


2001 April 17 


2017 


0.601 


22 


9083 b 


2001 October 1 


2183 


0.683 


165 


8619 b 


2001 November 27 


2240 


0.712 


-131 


9083 b 


2002 January 19-20 


2294 


0.738 


-82 


9337 b 


2002 July 4 


2460 


0.820 


69 


9420 b 


2002 December 16 


2664 


0.921 


-115 




2003 February 12-13 


2683 


0.930 


-57 




2003 March 29 


2727 


0.952 


-28 




2003 May 5 


2764 


0.970 


27 




2003 May 17 


2776 


0.976 


38 




2003 June 1 


2792 


0.984 


62 




2003 June 22 


2813 


0.995 


70 


9973 b 


2003 July 5 


2825 


1.001 


69 




2003 July 29 


2851 


1.013 


105 




2003 September 22 


2904 


1.040 


153 




2003 November 17 


2961 


1.068 


-142 




2004 March 6 


3071 


1.122 


-29 


Note. 


All HST STIS 


spectra used 


in this analysis 


arc available on 



http: / / archive . stsci . edu/prepds/etacar/ or http : // etacar.umn.edu/ 



" Phase relativ e to the beginning of the RXTE low-state 1997.9604: JD2.450, 799. 792+2024 X <j> 
JCorcoraril2005l) b HST GO Program - K. Davidson (Pl) c STIS GTO Program - T. R. Gull (PI) 

2. OBSERVATIONS 

Spectroscopic observations were conducted with the 
HST STIS beginning 19 98 January 1, jus t after the X-ray 
drop, as seen by RXTE (|Corcoranll2005( ), and continuing 
at selected intervals through 2004 March 6. The data 
set covers more than one full spectroscopic period (5.54 
years), with a focus on the 2003.5 event as a part of 
the ij Car HST Treasury project. The relatively short 
time span of the observation (6.2 years) compared to the 
spectroscopic period, restricts comparisons between long 
term trends and cyclic variations. To better understand 
the long term trends we need to have high spectral and 
spatial resolution spectra with good temporal sampling 
across the next event, occurring in early 2009. All ob- 
servations utilized the 52"x0'.'l long aperture in combi- 
nation with the G430M or G750M gratings, yielding a 
spectral resolving power of R~6000— 8000 and a signal- 
to-noise of 50—100 per spectral resolution element. As 
observations were desired at critical phases of the pe- 
riod, spacecraft solar power restrictions determined the 
slit position angle. Hence, the observations were done 
with a range of position angles. When possible, obser- 
vations were scheduled with slit orientation at position 
angle of —28° or 152° to include the spatially resolved 
Weigelt B and D blobs within the long aperture. With 
superb cooperation with the HST schedulers, this was 
accomplished several times throughout the 6.2 year ob- 
servational interval. 

The advantage of HST STIS over ground-based instru- 
ments is the high angular resolution, which permits ex- 
clusion of the bright nebular contributions of scattered 
starlight and nebular emission lines. While an ideal ex- 
traction of two pixels (O'.'lOl) would be desirable, reduc- 
tion issues due to small tilts of the dispersed spectrum 



with respect to the CCD rows can lead to an artificially 
induced modulation of the point-source spectrum. A 
tailored data reduction using interpolation was accom- 
plished with half-pixel sampling by Davidson, Ishibashi 
and Martin (http://etacar.umn.edu) minimizing the 
modulation. We determined that a six half-pixel extrac- 
tion (0'.'152) of the online reduced data provided an accu- 
rate measure of the stellar spectra while minimizing the 
nebular contamination. The 0'.'152 extraction covers a re- 
gion corresponding to ^350 AU at the distance of r\ Car, 
i.e. much larger than the major axis of the binary orbit 
(30 AU). The six half-pixel extraction excludes most of 
the radiation from the Weigelt blobs located 230—575 AU 
from r\ Car, characterized by narrow line emission. How- 
ever, even with the 0V152 extraction the wind lines are 
slightly influenced by some of the strong nebular emission 
features. The observed spectra used in this paper, with 
corresponding orbital phase 8 , are summarized in Table[TJ 
The complete HST STIS spectra for all phases are avail- 
ab le on [h ttp : / / archive . stsci . edu/prepds/et acar /| 
or |http : //et acar . umn . edu/] 

3. SPECTRAL ANALYSIS 

The spectrum of r\ Car contains a wealth of informa- 
tion. The spectrum is dominated by wind lines from 
rj Car A, in particular the H I Balmer and Paschen series 
and numerous lines from the singly ionized iron-group 
elements. A detailed discussio n of the optical spe ctrum, 
and its formation, is given by iHillier et all d200ll) . The 
UV w ind spectrum is described in IHillier et al.l (2001, 
2006), where the general characteristics are discussed in- 

8 All observations in this paper are referenced to the begin- 
ning of the RXT E low-state 1997.9604: JD2,450,799.792+2024x</> 
(Corcoran 2005) 
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Fig. 2.— He I A7067 profiles from 1998.0 to 2004.3. The He I 
lines profiles, mostly blueshifted relative to the —8 km s _1 system 
velocity, are asymmetric and vary in intensity and velocity with 
phase. The 2003 June 23 (Julian Day 2,452,813 = JD: 2813) ob- 
servation shows greater absorption at higher velocities (t>~— 800 
km s _1 ) than any of the other observations. This large blue-shift 
is due to a second absorption system at higher velocity. At phase, 
0=0.21, the absorption is weak and the velocity measurement is 
uncertain. 
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Fig. 4. — The Hel lines can be described by a broad P-Cygni 
wind profile with the superposition of several narrow emission com- 
ponents as exemplified with the He I A7067 profiles and component 
fits. Solid line: Observed spectrum. Dashed line: Fitted line pro- 
file based on a sum of the individual components (dotted lines). 
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Fig. 3. — Surface plot for Hel A7067. The line profile's variation 
with phase is presented with the intensity color-coded from black to 
red with increasing line emission. All spectra used in the figure has 
been normalized to the flux at —1200 km s^ 1 . The observed phase 
points are marked to the left in the figure. Notice that the two 
emission components merge just before periastron with a distinct 
velocity shift at phase, 0=1.0 (2003 June 29). The absorption, in 
contrast to that shown by Hi, is most prevalent at 0=0.8— 1.0. The 
contours represent the flux level at —1200 km s _1 plus 0.70 and 
0.85 of maximum line emission. 

eluding the absence of wind lines from r/ CarB. We con- 
centrate on the orbital phase dependance of four species 
whose spectra dominate the optical wavelength region: 
H I (ionization potential 13.6 eV), He I (IP 24.6 eV), 
Fe II (Fe II IP 16.2 eV with Fe III IP 30.7 eV), and [N II] 
A5756 (N I IP 14.5 eV with N II IP 29.6 eV). Due to 
the large abundance of hydrogen and its ionization po- 



JD: 3071 



ID: 2904 
0=1.04 

JD: 2851 
0=1.01 

JD: 2825 
0=1.00 

JD: 2813 
0=0.99 

JD: 2776 
0=0.98 

JD: 2460 
0=0.82 

JD: 2183 
0=0.68 

JD: 1623 
0=0.41 

JD: 1231 
0=0.21 

JD: 0891 
0=0.05 




-1000 -500 500 1000 

Velocity (km s" 1 ) 

Fig. 5. — Hi A4103 variability in intensity and velocity over the 
spectroscopic period. The emission centroid and emission fluxes 
are close to the system velocity (—8 km s ), and it and the emis- 
sion fluxes do not generally undergo the dramatic changes observed 
in He I. However, significant changes in the emission line radial ve- 
locity and the Hi emission line strength occurs near phase, 0=1.0. 
Changes are observed in the P-Cygni absorption over the entire cy- 
cle, with the absorption weakest at 0=0.21—0.98. The absorption 
velocity varies over the period but with a smaller amplitude than 
observed in He I. 
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Fig. 6. — Surface plot for Hi A4103. The line profile's variation 
with phase is presented with the intensity color-coded from black 
to red with increasing line emission. All spectra used in the figure 
has been normalized to the flux at —1200 km s —1 , and the phase 
points are marked on the left of the figure. The emission profiles is 
nearly constant, both in intensity and velocity, throughout most of 
the period. Significant changes are seen in the emission profile near 
phase, </>=1.0; there is a reduction in flux and a velocity shift. An 
abrupt strengthening of the P-Cygni absorption occurs after phase, 
0=1.0. The contours represent the flux level at —1200 km s _1 plus 
0.70 and 0.85 of maximum line emission. 

tential, H I lines form throughout the stellar wind struc- 
ture. Iron shows similar behavior to hydrogen but due 
to its ionization structure, Fe II lines form in the outer 
lower excited regions of rj Car A's wind, with the interior 
region dominated by more highly ionized iron. Helium is 
mainly neutral in the environment of the primary star. 
However, the energy level structure (see Figure [lj is such 
that high energy photons are required even to populate 
the lowest excited states. He I lines represent the highly 
excited regions of the stellar wind or the bow-shock. The 
[N II] A5756 line is a parity forbidden line, and is thus 
formed in the outer, low density, ionized region. 

STIS optical spectra of the star show that the excited 
He I stellar lines are peculiar relative to lower excitation 
wind lines: the He I lines a re blueshifte d relative to the 
system velocity (—8 km s _1 , [Smith 2004}) over most of the 
5.54 year cycle, unlike the lower excitation lines, whose 
emission is centered near system velocity. The He I P- 
Cygni lines in the 77 Car spectrum are characterized by 
a broad asymmetric P-Cygni profiles which often exhibit 
multiple, discrete peaks. Another unusual characteris- 
tic of the He I lines is the abrupt, large velocity shift 
towards the system velocity, exhibited in both emission 
and absorption, which occurred during the 2003.5 mini- 
mum. The variations of the He I P-Cygni lines in both 
strength and radial velocity suggest that most of the He I 
lines are influenced by t he ionizing flux of the companion 
star (jHillier et al.ll2006l ). 

While the spectrum of 77 Car is rich in wind lines, many 
lines are not suitable for detailed analysis. A large num- 
ber of He I lines are located in the STIS CCD spectral 
range (1600—10,300 A), but many, such as the commonly 
used He I A4473, are either blended or too faint for re- 
liable interpretation. Even with the 0'.'152 extraction, 
a few of the lines are influenced by emission from the 
surrounding nebula and the Wcigelt blobs. The He I 
AA5877, 6680 lines are slightly affected by narrow, neb- 
ular [Fe II] and [Ni II] emission lines, respectively. The 



least contaminated He I transition in the spectrum, for 
both absorption and emission, is He I A7067. Hence, 
we use this spectral line to illustrate the variability in 
absorption and emission. Unfortunately, the absorption 
in He I A7067 is often weak and difficult to measure. 
We have therefore also utilized the absorption portion of 
He I A5877. The He I A5877 line has the same lower state 
as He I A7067, but its oscillator strength is almost one 
order of magnitude greater and its absorption is corre- 
spondingly stronger. The He I A5877 emission is heavily 
affected by the strong Na I A5892 resonance absorption 
from the stellar wind, the intervening ejecta and the in- 
terstellar medium. The He I A5017 line is badly blended 
with Fe II A5020, and consequently little information can 
be gleaned about the emission profile. However, we can 
estimate the influence of the blend on the associated ab- 
sorption using other members of the same multiplet, such 
as Fe II A5170. We find that the absorption is dominated 
by the He I component throughout most of the period 
but is overwhelmed by Fe II absorption for a short in- 
terval after periastron (Section 14. ip . A complete set of 
lines used in this analysis is presented in Tabled while a 
detailed discussion of a subset of the least blended tran- 
sitions is presented in Section 4. 

To investigate the strongly phase dependent and com- 
plex profile variability we utilized a variety of techniques. 
To illustrate these we use He I A7067 which shows strong 
variability at all phases. At least four components can 
be identified in He I A7067: one absorption component, 
two narrow emission components which give rise to the 
bulk of the observed emission, and a much weaker, broad 
emission component. All these components show signifi- 
cant time variability. 

To illustrate the observed spectral variations, we 
stacked continuum-normalized line plots in velocity space 
(Figure [2|). We also used a color-coded intensity plot to 
illustrate the variation as a function of orbital phase. 
We created the color-coded plots by stacking the spec- 
tra according to phase, normalizing to the flux at —1200 
km s _1 relative to the vacuum rest wavelengths of the 
wind lines, where the continuum is seen to be unaf- 
fected by the wind line variations. The flux between 
the observations has been calculated by interpolation 
in phase to help identify systematic trends. The nor- 
malization procedure is used for the visualizations only, 
and utilized to emphasize the spectral variations with 
phase. The systematic change in radial velocity is ap- 
parent, as is a discontinuity in the radial velocity of 
the He I absorption near phase 1.0 (Figure [3]). Quan- 
tifying the observed spectral changes is difficult, since 
both the shape and strength of the line change. In ad- 
dition, reliable measurements of absolute line fluxes are 
difficult since the amount of extinction is declining with 
time. Since 1997, the brightness of the central star, as 
observed by HST in Ha, has increas ed by over a fac- 
tor three ([Martin fe Koppelman|[2004T ). During the min- 
imum, s mall variations (10—20%) in continuum flux ley - 
els occur ([Martin fe Koppelman1l2004HFeast et al.ll2001h . 
and the variations are not likely to be associated with a 
variation in dust extinction. 

The P-Cygni wind profile of a single massive star can 
be described by emission originating throughout the wind 
with absorption of the background emission in line-of- 
sight towards the observer. For a binary system the wind 
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profiles are more complex due to excitation and motion 
of the companion star. In the case of 77 Car, lines such as 
in H I and Fe II are formed throughout the dominating 
primary wind. In the absence of a companion star, the 
He I P-Cygni lines must form in a high-excitation zone 
close to the primary star. With a companion present, 
significant He I emission can arise in an ionized region 
near the bow-shock, and hence at much larger radii. The 
structure is even more complex due to the highly eccen- 
tric binary orbit (e~0.9). The Coriolis effects lead to a 
highly distorted emission structure, especially near peri- 
astron passage. By contrast the He I absorption arises 
from a confined but asymmetric region in line-of-sight 
from the primary stellar source. Consequently, we have 
chosen to use the P-Cygni absorption component in our 
analysis. To determine the radial velocity variation of 
the absorption components we have fitted multiple Gaus- 
sians to the line profiles, as illustrated in Figure [4] With 
this fitting procedure additional emission components are 
needed, which do not necessarily indicate spatially dis- 
tinct emitting regions. The specific number of the narrow 
components used in the fit are chosen to mimic the gen- 
eral shape of the emission line and, as far as possible, to 
decouple the emission from the absorption. The num- 
ber of narrow components needed to fit each line profile 
changes with phase, but is consistent for all lines at a 
given phase. The accuracy of the measured radial veloc- 
ity is phase dependent, due to a weak absorption or the 
possible presence of multiple components. Typical resid- 
uals to multi-Gaussian fits were less than two percent. 

4. LINE PROFILE VARIABILITY 

The line profiles in the r\ Car spectrum exhibit a wide 
variety of shapes and variability. Below we discuss the 
general characteristics and observed changes. 
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Fig. 7. — Top: Radial velocities measured for the He I absorp- 
tion (diamond/dashed) and Hi (crosses/solid). Note the dramatic 
velocity shift for the H e I lines after pe riastron (</>=1.0). Bottom: 
X-ray light-curve from l lCo rcoran 20Q5|) . 
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4.1. He I Profiles 

The He I line profiles are complex and appear to be a 
composite of multiple features. The profiles change con- 
siderably between observations, especially in the phase 
interval O.85<0<1.2O (Figures [2] and [3]). On top of a 
broad weak emission are relatively narrow emission com- 
ponents that change rapidly with time. All spectral line 
characteristics (i.e. strength, shape, and radial veloc- 
ity) vary systematically with phase. For each observa- 
tion, the He I profiles are shifted to the blue relative to 
the position of the H I lines (compare Figures [2] and [3] 
with Figures [5] and EJ. While both the H I and He I 
emission weaken during the minimum, the velocity shift 
for the He I emission is pronounced, and much larger 
than the corresponding shift observed in H I (see Sec- 
tion 14. 2p . We observe these shifts especially on the red 
side of the emission line where the profile is unperturbed 
by the blueshifted absorption. With phase, the He I 
emission increasingly shifts to the blue until just before 
the X-ray low-state begins. Between phases 0.995 and 
1.013, the intensity of the emission abruptly drops and 
the centroid of the emission line shifts to more positive 
velocities. After the minimum the emission is gradually 
more blueshifted. At some epochs the profile is nearly 
symmetric while at others the profile shows two or more 
narrow peaks. The limited sampling of this varying pro- 
file, and possibly spatially limited resolution, complicates 
the interpretation of the emission line radial velocities. 
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Fig. 8. — Radial velocities for the He I absorption. The typical 
measurement error is ~20 km s — 1 . Data points recorded between 
and 1.2 are re-plotted one period later to show periodic variability. 



Detailed modeling and possibly higher spatial/spectral 
resolution are required for a quantitative interpretation 
of the velocity variations. 

The He I absorption components vary in velocity and 
strength and is traceable over rj Car's entire 5.54 year 
spectroscopic period. The absorption velocity varies sys- 
tematically (except between </>=0. 0—0.2) with a discon- 
tinuous jump of more than 200 km s _1 near periastron 
(Figure Ej). 

In Figure [8] we show the measured absorption radial 
velocities for He I AA5877, 6680, 7067, 7283. The be- 
havior previously discussed can be observed in all the 
lines. However, there is scatter in the radial velocity 
measurements of different lines at the same phase. The 
absorption in He I and H I, apparent in Figures [3] and 
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Spectral Lines Used in This Analysis 
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Note. — Energy levels and wavelengths for Hi and the Hei triplets are weighted 
averages due to multiple unresolved transitions. The log gf is a sum over all transitions 
in the multiplet. All wavelengths arc in vacuum. 
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Fig. 9. — Top: Emission equivalent widths for Hei AA6680, 
7067. Middle: Absorption equivalent width for He I AA6680, 7067. 
Hei A6680 shown by asterisks and He I A7067 by diamonds. Bot- 
tom: Difference in absorption equivalent width as a function of 
phase. Error bars on the data points represent a 20% uncertainty 
in measurements for each of the two lines. The shift over the spec- 
troscopic cycle indicates a change in dominant excitation mecha- 
nism. During periastron an increase of equivalent width in Hei 
A7067 (triplet) relative to He I A6680 (singlet) indicates more ion- 
ization/recombination relative to photoexcitation. 

[6l is strongest at different phases. The He I absorption 
is strongest in the phase interval 0.8<</><1.0, and rela- 
tively weak at other phases. Conversely, the H I P-Cygni 
absorption, particularly of the Balmer series members, 
strengthens after 0=0.98, is extremely strong at 0=1.0, 
a nd weakens afte r the m inimum. 

Damincl i et al.l ((2000) investigated the variability in 
He I A6680 emission using ground-based, seeing-limited 
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Fig. 10.— Surface plot for Hi A3836. In contrast to Hi A4103, 
P-Cygni absorption is readily detectable at all phases. Like Hi 
A4103, and in contrast to the Hei lines, the Hi A3836 emission 
profile and radial velocity is relatively constant. Near 0=1.0, the 
emission profile weakens significantly with an associated velocity 
shift. The contours represent the flux level at —1200 km s _1 plus 
0.70 and 0.85 of maximum line emission. 

spectra. They concluded that the changes with phase 
were caused by excitation effects rather than motion 
of either of the objects in the binary system. The 
vel ocity amplitude measu red in He I A6680 emission 
by iDamineli et alJ (|2000t ) is in qualitative agreement 
with our observation s. However, the results from the 
IDamineli et al.l (|2000D analysis are based on variability 
in emission line strength while our results originate from 
absorption line radial velocities. The results may not be 
directly comparable but do show similar results regard- 
ing the amplitude and overall characteristic of the radial 
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velocity curve. 

The relative strength of the He I triplet lines to the 
singlet lines (see Figure provide information whether 
the lines are formed by photoexcitation or photoioniza- 
tion/recombination processes. An important question for 
understanding the origin of the He I absorption is what 
determines its strength: is it primarily due to optical 
depth effects or does the fractional coverage of the con- 
tinuum emitting source also play a role? 

We comment on the He I lines that can or cannot be 
studied in the r\ Car spectra covered by STIS: 

1. Absorption from the metastable state ls2s 3 P. All 
observable transitions from this term are severely 
blended. He I A3889 is, for example, blended with 
H I A3890. 

2. Absorption from ls2s 1 S, such as He I A5017. This 
line is blended with Fe II A5020. However, the ef- 
fect of the Fe II line can be determined using Fe II 
A5170 as Fe II AA5020, 5170 have the same lower 
energy state. Since the gf value for Fe II A5170 is 
similar to that of Fe II A5020, absence of signifi- 
cant absorption in Fe II A5170 for any specific ob- 
servation indicates that the absorption associated 
with Fe II A5020 is due to the He I line. From 
0=0.2—1.0, the He I contribution dominates the 
absorption, which is reasonable since the ls2s : S is 
a long-lived state. Late in this interval, the absorp- 
tion depth exceeds 50%, indicating that a signifi- 
cant fraction of the background emitting region is 
covered by He I gas responsible for the absorption. 
Unfortunately, from 0=0.0—0.2 and 1.0—1.2, the 
Fe II component becomes strong while, as shown by 
other He I lines, the He I A5017 absorption weak- 
ens. Hence, He I A5017 cannot be used to trace the 
He I absorption across the minimum. Absorption 
associated with He I is seen on He I A3966, which 
is a similar transition as He I A5017, but this line 
is severely blended. 

3. Absorption from ls2p 3 P, including He I AA4473, 
5877, 7067, where He I A5877 is the strongest. The 
strength of the absorption lines are partially con- 
trolled by optical depth effects, and not dominated 
by effects related to partial coverage of the contin- 
uum source. 

4. Absorption from ls2p 1 P, including He I AA4923, 
6680, 7283. He I A6680 is the least contaminated, 
but He I A7283 is usable in this analysis. 

Analysis of the He I absorption shows, that at some 
phases, the He I in the He + zone absorbs much of the 
flux from the back-illuminating source, r] Car A. For ex- 
ample, on 2002 July 4 (0=0.820), the central intensity in 
He I A5017 is less than 50% of the adjacent continuum 
level. We note that the line absorption equivalent widths 
do not vary uniformly. There is a significant difference 
in variation between lines as shown in Figure [9] where 
the equivalent widths of He I A6680 (singlet) and He I 
A7067 (triplet) are compared. A significant change in 
their relative strengths occurs around 0=0.9—1.1, indi- 
cating changes of the physical conditions of the absorb- 
ing regions. Figure[9]shows a greater equivalent width for 



He I A7067 near periastron compared to He I A6680, indi- 
cating ioniasation/recombination to be the major mech- 
anism for populating the He I energy states when rj Car B 
is close to i] Car A. 

We created radial velocity curves for the He I absorp- 
tion (Figure EJ) in He I AA5877, 6680, 7067, 7283, and a 
mean radial velocity curve (Figure [7]) . As with the peak 
absorptions, the most blueshifted of He I and H I occur 
at different phases. The greatest blueshift for He I occurs 
leading up to the X-ray low-state (0=1.0). In compari- 
son, the greatest blue-shifts for H I and Fe II occur during 
and post minimum and are noticeably less than observed 
in the He I lines. 

4.2. HI Profiles 

We use H I A4103 (US) to best illustrate spectral vari- 
ations in hydrogen. This is a moderately intense, un- 
blended line and is influenced by radiative transfer ef- 
fects in a less complicated way than Ha. The Hi5 line, 
as presented in Figures [5] and El shows phase dependent 
variations in both emission and absorption. The US emis- 
sion component varies substantially less than observed 
in He I (see Section 14. 1|) . Normalized to the contin- 
uum flux, the H I lines do not vary greatly in intensity, 
although a change occurs close to periastron (0=1.0), 
where the peak emission in H5 drops by approximately 
30%. This behavior has previously been discussed by 
iDavidson et al.1 (|2005h . where it was interpreted as ev- 
idence of a shell-ejection event. During the minimum, 
the US profile is flat topped (0=1.01 in Figured]), un- 
like profiles observed at other phases. One reason for 
invoking a shell-ejection is that a similar line profile was 
not observed during earlier spectroscopic minima. An 
alternative explanation is that 77 Car exhibits long term 
mass loss variations, and that these are superimposed 
on the phase related variations. The velocity variation 
of the emission is small (<25 km s _1 ), except between 
phases 0.9 and 1.1, and is illustrated by the stable be- 
havior of the red side of the emission profile. The ve- 
locity shift observed on the red portion of the H I lines 
is, however, dependent on a good understanding of the 
continuum and the spectral line shape. The measure- 
ments of the position of the absorption component yield 
a velocity ampl i tude o f ~60 km s — 1 . In an earlier study 
iDamineli etall (|2000fl . using ground-based spectra, de- 
rived a radial velocity curve using a large number of ve- 
locity measurements in H I Pa7 and Pa5. They derived a 
velocity amplitude of 50 km s , centered at the systemic 
ve locity. However , their orbital solution was, according 
to IDavidson et al.l (|2000D , inconsistent with higher spa- 
tial resolution HST spectra. A major difficulty in inter- 
preting the velocity shifts in the hydrogen lines occurs 
because the P-Cygni profile is confused in ground-based 
data, by nebular scattered radiation and emission from 
the Weigelt blobs. 

The major part of the H I emission is formed in the 
extended wind of 77 Car A. Some of the variations and 
structure seen on top of the broad emission may arise 
from asymmetries in the primary wind, or alternatively 
from changes in emission from the bow-shock associated 
with the wind-wind collision. Similarly, velocity varia- 
tions could be due to the motion of 77 Car A, or to changes 
in the wind structure. In both cases the observed ampli- 
tude of the velocity variations will be reduced since the 
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Fig. 11. — Illustration of Fen A5170 variability in intensity and 
velocity over the spectroscopic period. The absorption component, 
and emission on the blue side, are strongly variable with phase. 
Conversely, emission on the red side is relatively stable in strength, 
and in radial velocity. The contours represent the flux level at 
— 1200 km s — 1 plus 0.7 and 0.85 of maximum line emission. 

H I emitting region is extended well beyond r\ CarB's or- 
bit. Consequently, i t takes a finite time to replenish the 
gas. In the model of iHillier et ail (|2006[) . roughly 50% of 
the emission originates beyond i?^1500i?©. A character- 
istic flow time (R/va) is thus 25 days, assuming 1^=500 
km s _1 which is rj CarA's terminal wind velocity. Re- 
combination time scales are much shorter. 

In contrast to the emission, the H I P-Cygni absorp- 
tion is highly variable. At some phases (e.g. just after 
0=1.0) the absorption is strong, while at other phases 
(O.5O<0<O.95) it is much weaker, and in some lines dif- 
ficult to detect. For US the most striking behavior of 
the P-Cygni absorption is the rapid strengthening near 
0=1.0. At 0=0.984 (2003 June 1) no absorption can be 
seen in while it is marginally detectable on Ha. By 
0=0.995 (2003 June 22), the absorption increased dra- 
matically. P-Cygni absorption can be seen on high-n 
members of the Balmer series (Figure [T0|) at all epochs 
observed with HST STIS. 

The H I absorption shows radial velocity variations 
across the period (Figures [HI El an d fTU)) . During most 
of the period, the H I absorption velocity typically is 
between —330 to —400 km s" 1 , but varies substantially 
around periastron. In the phase interval O.95<0<1.2O, 
the absorption radial velocity falls from approximately 
—400 to —500 km s _1 , before recovering to its pre- 
minimum level. 

4.3. Fell Profiles 

The Fe II profiles also show intensity and radial veloc- 
ity variability. To illustrate the Fe II variability we use 
Fe II A5170 (see Figure [TT]) which is one of the stronger 
Fe II wind lines in the optical spectrum of r] Car. Ac- 
cording to models of luminous blue variables (LBVs) like 
•q Car, the Fe II emission is primarily produced b y con- 
tinuum fluorescence (e.g. IHillier et al.lll998ll2001h . This 
line provides a useful comparison to the He I A5017 + 
Fe II A5020 feature (Figure [12]), The third member of 
the multiplet, Fe II A4925, shows a similar emission pro- 
file as Fe II A5170 at all phases. 

Radial velocity variations are difficult to measure since 
the profile shows marked variability on the blue side. 
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Fig. 12. — Illustration of He I A5017 which is blended with the 
stronger Fell A5020. The velocity scale is for He I A5017. Fen 
A5020 has the same lower level as Fen A5170 (Figure [TTJ , and be- 
longs to the same multiplet. Comparison shows that absorption in 
the He I component dominates between ci=0.2— 1.0. The contours 
represent the flux level at —1200 km s — 1 plus 0.70 and 0.85 of 
maximum line emission. 

While some of this variability appears to be directly re- 
lated to the P-Cygni absorption, much of it can be as- 
sociated with change in the emission on the blue side. 
Conversely, the red side of the profile, with the exception 
of the spectrum recorded at phase 0.21 (1999 February 
21) shows minimal variability, with the velocity changes 
typically being less than 20 km s _1 . Similar to the H I 
line profiles, the Fe II absorption components are weak 
throughout most of the 2024 day cycle. They begin to 
strengthen around 0=0.98, are strong just after perias- 
tron (0=1.0), but have virtually disappeared at 0=1.2. 

Given the strength of the Fe II emission in the 77 Car 
spectrum, the weak P-Cygni absorption in all profiles 
of lines in the multiplet, including Fe II AA5020, 5170, 
throughout the broad maximum is difficult to explain 
with a spherically symmetric model. At times, the Fe II 
emission is double-peaked, while at other occasions the 
blue side of the profile is severely weakened relative to 
the red side. 

The spectral behavior observed in Fe II can indicate 
distinct emitting regions, one of which is highly variable. 
Alternatively, the bulk of the Fe II emission could come 
from the primary stellar wind. Variations would then be 
explained by significant perturbations to the ionization 
structure of the wind region that affects the blueshifted 
emission. Such a scenario is consistent with a hot com- 
panion star ionizing the wind on the observer's side of 
77 Car A. 

4.4. [N II] Profiles 

The [N II] A5756 line shows a complex structure with 
multiple narrow emission peaks spanning velocities of 
-500 to +500 km s^ 1 (Figures QJ& HU- This line profile 
is significantly broader than the He I emission lines and 
is centered about the —8 km s _1 system velocity. Near 
0=1.0 the emission is weak, with the observed line repre- 
senting intrinsic [N II] emission from the primary wind. 
At other phases the emission is stronger, more structured 
with the emission weighted towards the blue. Such emis- 
sion can arise due to additional ionization of the primary 
wind by rj CarB, or line formation in the bow-shock. At 
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Fig. 13. — [Nil] A5754 variations over the spectroscopic period. 
The line profile is asymmetric and can be described with four nar- 
row components which mutually vary in intensity but not in veloc- 
ity with phase. The line is centered near at the —8 km s -1 system 
velocity. There may be a broad underlying component associated 
with the stellar wind. 



phase 0.21 and 0.68 the [N II] profiles differ substantially 
from those in He I. The emission is more heavily weighted 
to the blue side. The [N II] emission components show 
little or no velocity variation; only the intensity of in- 
dividual components changes with orbital phase. The 
lack of velocity variation indicates that the [N II] emis- 
sion is formed beyond the inner, denser portions of the 
primary wind or the wind-wind collision region. The 
[N II] emission, as with numerous [Fe II] lines and other 
forbidden line species, originates in lower density, o uter 
wind regions. Using the model of Hil lier et al.l 1)200 If ) the 
extended [N II] emitting region coincides with, and ex- 
tends well beyond, the binary orbit (semi-major axis ~15 
AU). While the critical density of the 1 T) state is approx- 
imate ly 8.6 x 10 4 cm~ 3 at 10 4 K (jOsterbrock fe Ferlandl 
2006) , most of the emission may come from larger parti- 
cle densities if hydrogen has already recombined. 

5. SINGLE STAR SHELL-EJECTION VS. BINARY DRIVEN 
VARIABILITY 

A proper interpretation of the line profile variability 
depends on the excitation mechanisms in the line forma- 
tion region. Since rj Car is considered to be an LBV, a 
natural interpretation is that the observed emission and 
absorption line variations are evidence for periodic shell- 
ejections. Sh e ll-ejec tions were invoked, for example, by 
IZanella et alj (|1984f ) to explain the observed variation 



of high-excitation features seen in ground-based spectra. 
The discovery of the 2024 day period suggested that the 
variations could be explained by orbital modula tion in a 
binary system (e.g. iDamineli et al.l lT997. 199<|). In this 
model, flux from rj Car B ionizes more completely a por- 
tion of the outer wind on the side of rj Car B and weak- 
ens the hydrogen absorption. For this model to work, 
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Fig. 14. — Illustration of the [Nil] A5756 variability in intensity 
and velocity over the spectroscopic period. The contours represent 
the flux level at —1200 km s -1 plus 0.70 and 0.85 of maximum line 
emission. 

our line-of-sight must intersect the ionized zone 9 . At pe- 
riastron the flux from rj CarB is absorbed/occulted by 
rj Car A and its inner wind, and the H I absorption lines 
are stronger. 

A single star shell-ejection model is still invoked by 
some authors to ex plain the observed lin e and contin- 
uum v ariations fe.g. iDavidson et al.ll2005t ). ISmith et al.l 
(2003) invoked shell-ejection to explain the observed vari- 
ation of the H I P-Cygni profiles. In their model, the cen- 
tral source has a fast polar wind which is denser than the 
equatorial wind. At most phases the equatorial wind is 
highly ionized and the H I absorption is weak. However, 
during the spectroscopic event a shell is ejected and the 
wind becomes symmetric with enhanced H I absorption. 
The systematic variation of the H I absorption profiles 
seen in radiation reflected by dust in the Homunculus and 
the symmetry about the polar axis, support this model. 

Both models can potentially explain the observed vari- 
ations of the Weigelt blob spectra. In the single star 
shell-ejection model, radiation from the central source 
preferentially escapes in the equatorial regions, ionizing 
the Weigelt blobs. During the shell-ejection, this radia- 
tion is blocked, and hence the ionization of the Weigelt 
blobs decreases. In the binary model, rj CarB supplies 
the harder photons necessary to produce He I, Fe III etc. 
When rj Car B approaches periastron, rj Car A's extended 
wind and perhaps rj Car A itself, shields the blobs from 
the ionizing radiation. 

Can we distinguish between the single star shell- 
ejection model and the binary model? We interpret the 
present spectroscopic data as strong support for binarity 
over a single star shell-ejection: 

1 . The He I emission lines are consistently blueshifted 
throughout most of the 5.54 year period. This 
is not explainable with an axial-symmetric shell- 
ejection which would give equal amounts of 
blueshifted and redshifted emission. It is naturally 
explained by the binary model. In the binary model 
the He I emission is produced by ionization by the 

9 By ionized we mean H/H~'~~10 — 6 . In this case the Balmer P- 
Cygni absorption tends to be weak. When H/H + ~10 -3 or larger, 
the absorption is strong. 
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far-UV flux from r\ Car B followed by recombina- 
tion. The orbit is highly elliptical and at apastron 
(and for most of the orbit) rj CarB is in front of 
rj Car A, and ionizes the portion of 7/ Car A's wind 
flowing towards the observer. Therefore, for most 
of the orbit, the He I is formed in the blueshifted 
part of 77 Car A's wind. 

2. The double peak profile seen in the He I emis- 
sion (see Figure [2]), and the velocity variations of 
these peaks, are qualitatively consistent with emis- 
sion arising from spatiall y separated r egions near a 
colliding wind interface (|Luehrsl[l997f) . The shape 
and behavior of the He I emission lines in the 77 Car 
spectrum are similar to the line variations observed 
in other colliding wind binaries. The C III A5698 
line is observed to show dramatic line p rofile vari- 
ation in the spectra of W R 42, WR 48 (iHill et al. 
2000 L [2002h and WR 79 dLuehrsI H9971 kill et al. 
200dl2002t ). C III A5698 is not present in the 77 Car 
spectrum. 

3. As noted earlier, a striking characteristic of the He I 
emission is the dramatic velocity shift across peri- 
astron passage. This is similar to radial velocity 
variations of stars in highly eccentric binary sys- 
tems, and suggests the emission follows the orbit of 
one of the stars. In general it is difficult to see how 
a single star shell-ejection model can explain the 
complex behavior of the observed He I absorption 
profiles (see Section lf).2.1[) . Potentially, the velocity 
jump near 0=1.0 can be caused by a shell-ejection, 
with the absence of the higher velocity absorption 
from the previous shell disappearing because of re- 
combination. In a single star shell-ejection model, 
radial velocity variations in absorption would oc- 
cur since the ejected material accelerates to more 
negative velocities as it moves away from the star. 
However, in this model it is hard to explain why 
the most rapid decrease in velocity occurs close to 
2003.5 minimum rather than just after the shell- 
ejection, which would be necessary to provide sig- 
nificant absorption for the entire 5.54 year cycle. 
Furthermore, it is not easy to explain why the He I 
absorption increases in strength from 0=0.80 in a 
single star shell-ejection model, since this is well 
before a shell is ejected. 

4. While the increase in strength of the H I absorp- 
tion can be explained by a single star shell-ejection 
model, it is difficult to understand why the He I 
absorption has its greatest blue-shift at 0=1.0. In 
most single star shell-ejection models, the shell pro- 
duces low-velocity absorption just after ejection. 
The observed velocity change occurs across = 1.0 
and thereafter requires that the shell be ejected at 
a very high velocity (^—600 km s _1 ) and decel- 
erate just after 0=1.0 (to ^—300 km s^ 1 ) before 
re-accelerating. There is no obvious physical ex- 
planation for such behavior, other than the radial 
velocity variations are caused by orbital motion. 

5. Another concern is that the H I emission weakens 
near the time of the shell-ejection, before recov- 
ering to pre-2003.5 levels, but ejection of a shell 
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Fig. 15. — Radial velocity as a function of phase derived from 
the He I absorption. The solid curve shows a fit to the data with 
the parameter given in Table [3] (K= 140 km s _1 , e=0.90). Upper 
panel: Variation in the velocity amplitude, K=100 km s — 1 (dot- 
ted) and K=180 km s — 1 (dashed). Lower panel: Variation of the 
eccentricity, e=0.88 (dotted) and e=0.92 (dashed). The error bars 
for the measurements represent an estimated uncertainty of the fit 
and is proportional to the inverse of measured equivalent width. 

should produce more H I emission. In the binary 
model the weakening can be associated with the 
disruption of the primary wind by the companion 
star close to periastron. 

6. While high-n members of the Balmer lines show sig- 
nificant variations, P-Cygni absorption is present 
for these lines at all phases, implying that 77 CarB 
has a greater influence on the outer wind where, 
for example, Ha is formed. 

7. The variability in Fe II can be natura lly explained 
by a binary model. In the model of iHillier et al.l 
(2006), Fe II emission originates at large radii, ex- 
tending beyond the binary orbit (30 AU). The com- 
panion's far-UV emission will enhance the ioniza- 
tion of nearby outflowing material, weakening both 
the emission and absorption. Because the orbit is 
highly elliptical, for most of the orbit 77 Car B is on 
the observer's side of 77 Car A. On the other hand, a 
single star shell-ejection model can explain some of 
the observed variations if the ejected shell obstructs 
UV photons from 77 Car A, lowering the degree of 
ionization and enhancing the strength of the Fe II 
P-Cygni absorption. However, an axial-symmetric 
shell-ejection usually causes emission variations on 
both the blue and red sides of the line profile. Only 
strong variability on the blue side is observed. 

There are still uncertainties and some of the individual 
points raised above are debatable. However, the consis- 
tency between the different line profile diagnostics, the 
qualitative variation of the observed X-ray fluxes, the 
variation of the line intensities in the Weigelt blobs, and 
the location of the Weigelt blobs on the same side of 
r\ Car A as 77 CarB at apastron, strongly support the bi- 
nary scenario. Furthermore, the data suggest that the 
observed variations are a direct consequence of the vari- 
ability of the UV radiation and wind of r\ Car B. 

6. INTERPRETATION OF THE VELOCITY VARIABILITY 
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We have established that line profile variations ob- 
served over rj Car's spectroscopic cycle are consistent 
with a binary model. However, a quantitative analysis of 
the observed variations is difficult because the lower exci- 
tation lines arc formed in an extended part of the system. 
An exception is the He I absorption, that presumably is 
formed in a localized region close to the wind-wind inter- 
face. The radial velocity variations in the He I absorp- 
tion appear to be reminiscent of radial velocity variations 
from a star in a highly eccentric orbit: a gradual change 
in radial velocity through most of the cycle, with a rapid 
shift near periastron passage. We therefore explore what 
these velocity variations tell us about n Car A, rj CarB, 
and the interacting gas. 

We show how the He I absorption line variations can 
be described by a standard eccentric binary radial veloc- 
ity curve, and demonstrate that the derived value of the 
eccentricity is consistent with the value derived from the 
X-ray light-curve. To interpret the absorption compo- 
nents and the meaning of the velocity curve we present a 
detailed discussion of the colliding wind model and dis- 
cuss the origin of the He I emission and absorption com- 
ponents. 

6.1. The He I Radial velocity curve 

Line profiles recorded at several phases indicate that 
the He I absorption is a composite of multiple compo- 
nents. However, only the strongest, least blueshifted fea- 
ture, is traceable over the entire spectroscopic cycle. We 
created a radial velocity curve for the He I absorption, 
based on this component. A second high velocity absorp- 
tion component appears near the spectroscopic minimum 
(see Figure [4 ]) and disappears during the recovery phase. 
iStahl et al.l (|2005h noted a higher velocity component 
that appeared near the time of the X-ray low-state, and 
claimed it was a result of a shell-ejection. Over the brief 
period in which it was observable it followed a velocity 
pattern similar to the less blueshifted component, sug- 
gesting formation in adjacent regions. We have compared 
the velocity variation with the ra dial velocity curves for 
He I and H I emission derive d bvlDamineli et all (2000) 
and to the X-ray light-curve (|Corcoranll2005f) . Measured 
absorption velocities for the He I and H I are presented 
in Figure [7] along with a plot of the RXTE X-ray flux 
(2—10 keV). The maximum velocity of the He I absorp- 
tion is about —600 km s" 1 , which decreased to about 
—300 km s -1 on 2003 June 22, just before the spectro- 
scopic minimum. The striking difference between the 
velocity variation for the He I and H I wind lines is the 
drop in velocity after periastron ( 0=0.05) which follows 
the X-ray drop seen by Chandra (|Corcor an 2006). This 
discontinuity may be a result of absorption in the gas be- 
tween the primary star and the bow-shock, while the H I 
absorption originates from a volume extended beyond 
rj CarB's orbit. The variation of the He I absorption 
component is similar to that of a spectroscopic binary. 
Figure [T5l shows the observed He I absorption velocities, 
along with sample velocity curves for a range of eccen- 
tricities and velocity amplitudes. Note that the spectral 
lines in this analysis have absorption components which 
are blended with their emission components, making the 
line profile fitting difficult and decreasing the measure- 
ment accuracy. The radial velocity curves for He I, plot- 
ted in Figure [T5l were calculated using the standard for- 




FlG. 16. — A cartoon representing the standard view of the 
colliding wind system in rj Car, in which the shock-cone is wrapped 
around rj CarB (B). The more massive star rj CarA (A) has the 
slower but denser wind. The arrows show the flow of the gas along 
the bow-shock. The position of HST is appropriate for u;=270° 
and the system is likely tilted about 41° into the sky-plane. 



TABLE 3 
Radial Velocity Curve Solution 



Parameter 



Value 



p 


2024.0 days 1 


i 


41 o2 


7 


-430 km s" 1 


K 


140 km s _1 


e 


0.9 




270° 



1 ICorcoranl II2005T) 2 IDavidson et al.l <2001f) 

malism (|AitkeiJll96l ): 



dz 
~dt 



K(e cos w + cos (y + to)) + 7 



(1) 



where dz/dt is the radial velocity along the line-of-sight, 
e is the eccentricity, uj the longitude of periastron, v the 
true anomaly 10 and 7 the average velocity. K is: 



K 



2na sin i 



(2) 



where i the inclination, a the major axis of the orbit 
and P is the spectroscopic period with an assumed be- 
ginning at the X-ray low-state coinciding with periastron 
passage. All other parameters are initially free, although 
we constrai n e>0.8 based on th e solution of the X-ray 
light-curve (|Corcoran et al-l feOOl). In Figure [15] solutions 
obtained by varying e and K are plotted. The adopted 
parameters are given in Table [3] 

6.2. Orbital Dynamics 

In most published binary models of the rj Car system, 
•q Car A is defined as the primary star and is dominating 
both the mass and luminosity of the system, n CarB, 
is usually assumed to be the hotter, less massive star 
with a weaker, although much faster wind. Support 
for such a scenario comes from the analysis of th e X- 
ray (|Corcoranl[2006h. the farJJV jlping et al.ll2005f ) and 
UV/optical (jHillier et al .||200lLl2006D spectra in additio n 
to the modeling of the Weigelt blobs (|Verner et al.f 2005') . 

10 The relation between the true anomaly, v, and the eccentric 
anomaly, E, is given by tan (^) = J ji| tan (^). E is related to 



the phase as 2w4 



- e sin (E). 
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Though there is no direct evidence regarding the orbital 
inclination or the relation of the orbit to the geometry of 
the Homunculus, the usual assumption is that the incli- 
nation of the orbit is the same as for the symmetry axis 
of the Homunculus, and that the orbital plane lies in the 
intervening dis k between the l obes o f the Homunculus . 
Most models (|Corcoran et all 120011: ISmith et all l2004h 
assume that the longitude of periastron is such that, prior 
to periastron passage, 77 CarB is moving away from us 
and 77 Car A is moving towards us, as shown in Figure \W[ 
iFalceta-Goncalves et al.l (|2005l ) suggested that rj Car B is 
moving towards us just prior to periastron passage. It is 
more certain that the system has a large eccentricity and 
that, because of the short duration of the pan-chromatic 
variations, periastron passage occurs in the time interval 
of the X-ray/ionization low-state. 

The He I P-Cygni absorption components are the best 
localized features in the system since they require high 
energy photons and/or collisions to be produced and they 
arise in out-flowing material distributed along the line-of- 
sight to the more luminous star. Nevertheless, the inter- 
pretation of these features is not straightforward since 
it is not well understood where the absorbing material 
is located. The absorber is either associated with the 
wind of r\ Car A, with the colliding wind interface, the 
wind of 77 Car B, or a combination of these. Similarly, 
the He I emission can either arise in the bow-shock or 
in the portion of the primary wind ionized by 77 CarB. 
Downstream absorption in the bow-show has previously 
been observed in other binary systems such as V444 Cyg 
(jShore fc Brownl I1988D . However, the thickness of the 
bow-shock in the 77 Car system and the projection in line- 
of-sight provide little absorbing material between 77 Car A 
and the observer throughout the orbital period. The ion- 
izing photons of 77 Car B penetrate deeply into 77 Car A's 
wind creating a substantial He + zone. In the following 
sections we describe the possible location of the absorb- 
ing system, and discuss implications on the derivation of 
the system parameters. 

6.2.1. Velocity variations caused by orbital motion 

The velocity curve shown in Figure[T5]is very similar to 
a radial velocity curve for an object in a highly eccentric 
orbit. If Figure fT5l describes the orbital motion of one of 
the two objects in the 77 Car system then the blueshifted 
He I absorption over most of the spectroscopic period and 
the geometry of the system (including the position and 
variability of the Weigelt blobs), imply that we are trac- 
ing absorption in the dense primary wind. The terminal 
velocity of the absorption profile is consistent with the 
observed terminal velocity of 77 Car A and much lesser 
than the terminal v elocity of 77 CarB (^3000 km s _1 ) 
(jPittard fc Corcoran|[2003l ). Wind profiles with a termi- 
nal velocity as large as 3000 km s _1 are not observed in 
the 77 Car spectrum. 

The derived radial velocity curve (Figure IT5]) confirms 
the large eccentricity derived from the X-ray light-curve 
([Corcoran et al.ll2001h . but yields a surprisingly large K- 
value, K = 140 km s _1 . This value is inconsistent with 
the radial velocity amplitude measured in H I, whose ori- 
gin is associated with 77 Car A. To determine a mass of 
both stars in the binary system, radial velocity curves 
for both objects are needed. However, if only one of the 
objects is traced then a mass function describing the rela- 



tion between the two objects can be created. If the large 
He I velocity shifts during periastron are solely caused 
by the orbital motion of one of the stars, then the mass 
function is: 



/(mi,m 2 ) = 



K(P{1 - e 2 ) 



2\3/2 



2nG 



mi sin i 



m'A sin 1 



(mi + m 2 ) 2 



50M P 







using the values from Tabled and q=mi/m2 with the 
subscripts denoting the observed star (1) and the unseen 
star (2). With the exception of the possible detec tion of 
77 CarB in the FUSE spectrum (jlping et al.ll2005l ). there 
is little evidence for spectral features from the companion 
star. We therefore, initially, associate the He I absorp- 
tion with the primary star. If the observed velocities 
represent the motion of 77 Car A, a mass for 77 Car B for 
chosen values of 77 Car A's mass can be derived. This 
analysis is consistent with a mass of 77 CarB of 210Af Q 
for a ~20Mq mass of 77 Car A, i.e. 77 Car A is the lower 
mass object in the system with an extreme mass ratio, 
<7^0.1. This is the solution with the lowest total mass 
of the system, but is only one of many possible solutions 
with the set of parameters given in Table [3] The orbital 
plane is assumed to be in the same plane as the inter- 
vening disk between the Homunculus, z=41°, but this 
value is poorly determined. If using a limiting value of 
the inclination closer to 90° a different set of solutions 
with a higher mass ratio (0.4—0.6) is obtained, however, 
77 Car B still is the more massive object in the system. 
This is different than the conventional view of the system 
in which 77 Car A has a mass of ^1OOM and is the more 
massive object. 

If the conventional view, as we believe, is correct, it re- 
quires that the mass function we derive above is a gross 
overestimate of the true mass function. This could be 
due to underestimation of the eccentricity or inclination, 
or an overestimate of the velocity amplitude, K, or a 
combination thereof. Note, that the argument that the 
more luminous star is the more massive one, assumes the 
two objects to be unevolved. If this is not the case and 
for example 77 Car A is on the helium burning main se- 
quence, then 77 Car B may be the more massive object in 
the system. 

In the next section we describe how ionization effects, 
caused by 77 CarB, may influence the derived value of K. 

6.2.2. The Influence of Ionization on the Velocity Amplitude 

If the He I absorption is in states populated by recom- 
bination in a region of 77 Car A's wind, which is ionized by 
77 Car B and localized towards 77 Car A, then the observed 
absorption line velocities are a combination of 77 Car A's 
wind velocity in the absorbing region in addition to the 
velocity due to 77 Car A's orbital motion. Variations in 
the observed line velocities can either be produced by 
changes in the orbital velocity of the absorbing gas, in 
the characteristic outflow velocity of the gas, or a com- 
bination thereof. 

The location of the ionized region in 77 Car A's wind 
is difficult to determine accurately since it depends on 
the density profile through its distorted wind, combined 
with the size of the wind- wind collision shock- front which 
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(a) 





(b) 



Fig. 17. — Visualization of the He + zone located in the rj Car system. r\ Car A, the l arger red stella r component, is shown with a ~9 
AU diameter disk, consistent with the apparent 2 /an continuum diameter measured by Weigelt (2006). rj CarB, the source of the UV 
radiation, is the smaller blue disk. The gray structure represents the He + zone, the source of the He I emission and absorption. Top row 
(a) views the system with the plane lying in the paper: the Z-axis (red) points out of the paper, the X-axis (dashed arrow) and the Y-axis 
lie in the paper. The semi-major axis of rj CarB is in the Y direction and the semi-minor axis in the X direction. Bottom row (b) views 
the system with the Z-axis (red) pointed along the a xis of the bipola r Hom unculus with Z at —41°, rotated about —45° in position angle 
from north, consistent with the geom etry derived by [ Davidson ct al. (2001). Given the locations of the Weigelt blobs plus the orientation 
of the radio continuum monitored by Duncan & White (2003), the orbit is rotated into the disk plane such that the semi-major axis (Y 
direction) projects almost on top of the Z-axis. NOTE: this visualization does not take into account distortion due to the Coriolis force 
which is considerable near periastron (0=0.00), but small near apastron (</>=0.50). 

allows more of rj CarB's ionizing flux to penetrate into 
rj CarA's wind. We illustrate an undistorted He + zone 
for different phases throughout the 5.54 year period in 
Figure[T71 We note that even the two-dimensional model 
of the colliding wind region becomes highly distorted by 
the Coriolis force near periastron (J. Pittard, 2006, pri- 
vate communication) . Distortion of the wind- wind inter- 
face will greatly affect the derives radial velocity mea- 
surements, and the derived K value. Building a detailed 
three-dimensional model is beyond the scope of this pa- 
per. The model shown in Figure [T71 simply approximates 
the wind-wind interface and the He + boundary as two 
undistorted parabolical surfaces centered on the line be- 
tween the two stars. 

Figure [17] (top row) shows the system as viewed look- 
ing down on the orbital plane. We see that, other than 
near periastron (<t> — 0), the He + structure is predomi- 
nantly offset from rj Car A towards rj CarB. In this ori- 
entation, the velocity seen in emission will be centered at 
the system velocity, but the observer would see a range 
in velocities dependent upon the angular curvature of the 
bow-shock in three dimensions. The major part of the 
H I Lyman radiation escapes down the direction of the 
major axis and is directed in a small cone along the or- 
bital plane toward the Weigelt blobs. 

The bottom row in Figure [T71 shows the binary system 
viewed by the observer, assuming an inclination, i=— 41°, 
which is appropriate if the orbital plan is in the disk plane 
between the bipolar lobes of the Homunculus. The ori- 
entation depicted in the bottom row of Figure [17] is such 
that for most of the period, a significant portion of the 
He + zone overlaps the line-of-sight towards rj Car A, but 
at different radial velocities if the wind is accelerating 
outwards from rj Car A. 

Figure [18] qualitatively shows the effects of ioniza- 
tion on the material in front of the disk of rj Car A. 
The rectangular region represents the material towards 
rj Car A, and the dark grey region represents the He+ 
zone and scales with the He I column density, propor- 
tional to N(Hc + )xN(c~). Prior to periastron passage, 
the companion star mostly ionizes high-velocity wind 
material far from the primary in the direction of the 



Fig. 18. — A qualitative illustration showing the orbit of 
rj CarB in the reference frame of rj Car A and how the ioniza- 
tion of rj CarA's wind by r\ CarB can produce absorption line 
velocity variations. The dashed lines represent wind isovelocity 
contours. The dark grey area represents the section of rj CarA's 
wind which is ionized by rj CarB's flux. Prior to periastron pas- 
sage (top), rj CarB primarily ionizes high-velocity wind material 
far from r\ Car A. After periastron passage (bottom) rj CarB mostly 
ionizes low velocity wind material due to a buildup of cool wind 
material between the observer and rj CarB. This is a qualitative 
representation and not drawn to scale. 

shock-cone, which points away from rj Car A. After peri- 
astron passage, the shock-cone is tilted towards 77 Car A, 
so that r] Car B mostly ionizes low velocity material in 
rj Car A's inner wind. X-ray spectra also show a buildup 
of cool, absorbing material between the colliding winds 
and t he observer for a brief time after periastron pas- 
sage (jViotti et al.ll2002t ICorcoranl l2005f) . which blocks 
ionizing radiation from the companion star from reach- 
ing the outer, high velocity wind material. Thus, after 
periastron passage the observed velocity of the absorp- 
tion system is expected to shift to lower velocities. This 
means that ionization effects can produce larger velocity 
variations than would be observed from orbital motion 
alone, and therefore the velocity amplitude derived in 
Section 16.2.11 does not represent the actual velocity vari- 
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ation due to the orbital motion of rj Car A. If the ionized 
structure was known, the He I radial velocity curve could 
be corrected for the influence of ionization and the true 
orbital motion of 7/ Car A could be derived. However, 
this requires detailed modeling of the dynamics of the 
wind of 1] Car A and rj CarB. 

7. DISCUSSION 

Evidence from numerous observations suggest that 
•q Car A is orbited by a hot companion star in a highly 
eccentric orbit. The velocity curve derived from the He I 
P-Cygni absorption, and to a lesser extent, the variation 
of H I absorption show periodic behavior very similar to 
the radial velocity of a star in a highly eccentric orbit 
with a longitude of periastron w~270°, which is consis- 
tent with the orbital paramet ers derived from analy sis of 
the RXTE X-ray light-curve ([Corcoran et al.ll200ll ). and 
also confirm the variation seen in ground -based observa- 
tion of the Pa7 and Ya5 emission lines ([Damineli et al.l 
I2000T ) . Significant deviations from u> by more than about 
10° could not describe the observed radial velocity curve. 
Our solution req uires a much larger ecc entricity than the 
solution given bv lDamineli et alJ " i|200Ct ). e=0.9 instead of 
e=0.65— 0.85. This larger eccentricity helps resolve the 
discrepancy between the ground-based Paschen line ve- 
locities and an earlier report on P aschen line velocities i n 
two observations with HST STIS ([Davidson et al.ll2000f ) . 
In any event the more complete STIS dataset analyzed 
in this paper clearly show radial velocity variations of 
the type expected from orbital motion. Our analysis of 
the radial velocity curve shows that the orbit is oriented 
so that the semi-major axis is pointed toward the ob- 
server (oj=270°) but with some orbital inclination (likely 
i~41°) and that r\ Car A is in front of its companion at 
periastron. In our simplistic model the He I absorption 
forms in the wind of rj Car A due to ionization by the 
far-UV flux from rj CarB, and that the He + portion of 
rj CarA's wind approaches the observer prior to peri- 
astron passage, consistent with the increasing blueshift 
seen in the He I and H I lines. 

One surprising result of our He I absorption analysis, is 
the large derived value of the velocity amplitude, if =140 
km s _1 . Since the star which is in front at periastron is 
rj Car A, the large if- value along with the increasingly 
blueshiftcd absorption prior to periastron, would imply 
that rj Car A itself is the lower mass object, or, more pre- 
cisely, that the star whose wind dominates the system is 
the lower mass object. Such a scenario is not unusual. In 
most well-studied colliding wind binaries the lower mass 
star has the stronger wind. A case in point is WR 140, 
a canonical WC+O colliding wind system, which, like 
rj Car, has a long p eriod (P=7.9 years), h igh eccentric- 
ity fe=0.881±0.005 lMarchenko et all [20031 ). hard X-ray 
emission and in which the weaker-wind O star dominates 
the mass and luminosity of the system. An exception to 
this generalization is the LBV+WR system HD 5980, in 
which the more massive star is an eruptive variable that 
currently possesses the stronger wind though prior to its 
erupti on in 1994 the more mass ive star had the weaker 
wind ([Koenigsberger et al.ll2006l ). 

A more plausible explanation is that the if-value de- 
rived from the He I radial velocity curve is overestimated 
because of ionization effects, as discussed in Section l6.2.2l 
An exact modeling of the effects of rj Car B's far-UV flux 



on the ionization state of rj Car A's wind depends sensi- 
tively on knowledge of rj CarA's wind structure (which 
is distorted by the wind-wind interaction) and on the 
companion's far-UV flux, neither of which is well con- 
strained. We qualitatively represent these effects, but 
our model of an undisturbed primary wind being ionized 
by the far-UV of the secondary is too simple. Numer- 
ous effects may come into play. For example, the wind 
from rj Car A piles up on the colliding wind shock, while 
radiation and thermal energy from the wind shock may 
also excite helium. As the system approaches periastron, 
the evidence is strong that the secondary star approaches 
to within a few AU of rj Car A. When this happens the 
winds and atmospheres of the two stars may experience 
a major disruption. This could have a major influence on 
the observed line profiles and measured velocities. The 
Coriolis effect in this very eccentric orbit greatly influ- 
ences the structure and therefore, the derived if- value. 
It is a matter of controversy whether the colliding wind 
interface remains intact at this point. 

The He I and H I line profiles reveal important infor- 
mation about the binary system. However, only two ob- 
servations, the first about one week before the minimum 
and the second one week afterward, measure the rapid 
changes in velocity in He I and H I. Much more frequent 
observations with high spatial resolution are needed. We 
caution observers that because of nebular scatter of the 
stellar component and contamination of nebular emission 
lines, sub-arcsecond angular resolution is necessary. 

8. SUMMARY 

We have monitored the He I and H I emission and ab- 
sorption between 1998.0 and 2004.3. The He I lines gen- 
erally show blueshifted peculiar profiles, with a dramatic 
velocity shift across the spectroscopic minimum. These 
lines are likely produced by the flux of the hot compan- 
ion star, though energy released in the wind-wind inter- 
face region may also play a role. The He I lines show 
a multi-component structure, consisting of a broad P- 
Cygni feature with narrow peaks superimposed on top. 
The lines are generally blueshifted for most of the spec- 
troscopic/orbital cycle, showing that they are produced 
in the portion of the wind of rj Car A which is flowing 
towards the observer. We have interpreted the narrow 
emission line components as formed in spatially sepa- 
rated regions in the bow-shock. Hence, the He I lines 
forming in rj CarA's wind in the vicinity of the bow- 
shock, are the result of the far-UV flux from the hot 
rj CarB. The H I lines display homogeneous profiles cen- 
tered at the system velocity and show a smaller shift in 
velocity across the spectroscopic minimum. The smooth 
velocity curve for the H I lines implies line formation 
from a large part of the system. 

We have measured the radial velocity and equivalent 
width variations in the He I absorption and emission lines 
and created a radial velocity curve from the absorption 
lines. Using this radial velocity curve, we have calculated 
a mass function for the binary system and estimated the 
mass of the hidden companion star, and show that either 
the companion star is the more massive object in the sys- 
tem, or that our estimate of the velocity amplitude is too 
large due to ionization effects. We believe that the lat- 
ter explanation is the more likely one, and that detailed 
modeling of the radiative transfer in the system may al- 
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low us to constrain the mass function more accurately 
from the HST STIS data. 



The observations were accomplished through STIS 
GTO, HST GO and HST -q Car Treasury Team observa- 
tions and were made with the NASA/ESA HST Science 
Institute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under NASA 



contract NAS 5-26555. We thank the STScI schedulers 
who worked closely with T. R. Gull in designing the ob- 
servational programs. All planned observations were ac- 
complished at the requested intervals and virtually all 
observations were optimally exposed. All analysis was 
done using STIS IDT software tools on reduced data by 
K. Davidson, K. Ishibashi and J. Martin made available 
through the HST rj Car Treasury public archive. 
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